DESY 04-001 
IFT-0 1/2004 
ZU-TH-1/04 
^-ph/040104l| 



O : Three-loop matching of the dipole operators 

00 



X 



for 6 ^ §7 and b ^ sg 



Mikolaj Misiak^'^ and Matthias Steinhauser^ 



^ Institute of Theoretical Physics, Warsaw University, 

O : Hoza 69, PL- 00-681 Warsaw, Poland. 

O 

^ //. Institut fiir Theoretische Physik, Universitdt Hamburg, 
p i' Luruper Chaussee 1^9, D-22761, Hamburg, Germany. 

O , Institut fiir Theoretische Physik, Universitdt Ziirich, 

Winterthurerstrasse 190, CH-8057, Ziirich, Switzerland. 



Abstract 



We evaluate the three-loop matching conditions for the dimension-five operators that are 
relevant for the 6 — 57 decay. Our calculation completes the first out of three steps (matching, 
mixing and matrix elements) that are necessary for finding the next-to-next-to-leading QCD 
corrections to this process. All such corrections must be calculated in view of the ongoing 
accurate measurements of the B Xs7 branching ratio. 



1 Introduction 



The inclusive weak radiative i?-meson decay is known to be a sensitive probe of new physics. 
Its branching ratio has been measured by CLEO [1], ALEPH [2], BELLE [3] and BABAR [4]. 
The experimental world average [5] 

BR[B X,7, {E^ > j^ruf,)] = (3.34 ± 0.38) x 10"^ (1.1) 

agrees with the Standard Model (SM) predictions [6, 7] 

BR[B Xs-f, (E^>1.6GeV)] = (3.57 ± 0.30) x 10-^ (1.2) 
BR[B X,7, {E^ > ^mt,)] - 3.70 x 10"^ (1.3) 

Such a good agreement implies constraints on a variety of extensions of the SM, including 
the Minimal Supersymmetric Standard Model with superpartner masses ranging up to several 
hundreds GeV. These constraints are expected to be crucial for identification of possible new 
physics signals at the Tevatron, LHC and other high-energy colliders. However, any future 
increase of their power depends on whether the theoretical calculations manage to follow the 
improving accuracy of the experimental determinations of BR[i? — Xg'y]. 

As pointed out more than two years ago [6], the main theoretical uncertainty in the SM 
prediction for BR[i? Xg'j] originates from the perturbative calculation of the 6 — > 57 ampli- 
tude. It is manifest when one considers the charm-quark mass renormalization ambiguity [6] 
in the two-loop, Next-to-Leading Order (NLO) QCD corrections to this amplitude [7,8]. The 
only method for removing this ambiguity is calculating the three-loop, Next-to-Next-to-Leading 
Order (NNLO) QCD corrections. A sample NNLO diagram is shown in Fig. d 




Figure 1: One of the C(IO^) three-loop diagrams that we have calculated. 

Since rrify <^ M^r, such diagrams are most conveniently calculated using an effective field 
theory language. The electroweak-scale contributions are encoded into the matching conditions 
for the Wilson coefficients, while the fe-quark-scale contributions are seen as matrix elements 
of several fiavour-changing operators. Large logarithms ln(M^/mf) are resummed using the 
effective theory Renormalization Group Equations (RGE) that result from the operator mixing 
under renormalization. 
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The matching conditions and the matrix elements yield a renormalization-scheme indepen- 
dent contribution to the amplitude only after combining them together. Thus, both of them 
need to be evaluated to the NNLO. It is impossible to remove the charm-quark mass am- 
biguity by calculating the matrix elements only, even though the matching conditions are 
rric- independent . 

In this paper, we present our calculation of the three-loop matching conditions for the dipole 
operators {sLcr'^'^bjijF^i, and (si(T'"^T"6/j;)G'^^. All the other matching conditions that are rel- 
evant for 6 — > S7 at the NNLO originate from two-loop diagrams only, and were calculated 
several years ago [9]. Thus, our work completes the first (matching) step of the full C(a^) 
analysis of the considered process. 

The long history of the lower-order {0{1) and 0{as)) analyses has been summarized in 
Ref. [10]. As far as the NNLO calculations are concerned, fermion-loop contributions to the 
three- loop matrix element of the current- current operator {sc)Y_A{cb)Y_A are already known 
[11]. Three- loop anomalous dimensions of all the four-quark operators will soon be published 
[12]. Work at the remaining anomalous dimensions and matrix elements is in progress. 

In our present three-loop matching computation, we follow the procedure outlined in Ref. [9] . 
All the necessary diagrams are evaluated off-shell, after expanding them in external momenta. 
The spurious infrared divergences generated by the expansion are regulated dimensionally. 
They cancel out in the matching equation, i.e. in the difference between the full SM and the 
effective theory off-shell amplitudes. 

The scalar three-loop integrals are evaluated with the help of the package MATAD [13] designed 
for calculating vacuum diagrams. The fact that MATAD can deal with a single non-vanishing 
mass only is not an obstacle against taking into account the actually different masses of the 
W boson and the top quark. Expansions starting from = M\y and ^ Myy allow us to 
accurately determine the three- loop matching conditions for the physical values of and M^y. 

Our paper is organized as follows. In Section |21 we introduce the effective theory and give 
all the necessary renormalization constants. In Section 01 the unrenormalized one-, two- and 
three-loop SM amplitudes for b ^ s'j and b sg are presented. Section |3] is devoted to a 
discussion of the SM counterterms. The matching procedure is described in Sectional Explicit 
expressions for the resulting Wilson coefficients are given in Sectional We conclude in Section [7| 
Appendix A contains exact expressions for the coefficients of the expansions in (1 — M^/m'f) 
and M^/m'f. 



2 The effective theory 

Since our approach follows Ref. [9] very closely, we shall not repeat all the details given there. 
While the present article is self-contained as far as the notation is concerned. Sections 2 and 5 
of Ref. [9] are referred to for pedagogical explanations. 

The effective theory that we consider arises from the SM after decoupling of the heavy 
electroweak bosons and the top quark. Its off-shell Lagrangian reads 

Cefi = /:QCDxQED(n, d, s,c,b) + ^Y. IVclVcbC^ + v*v,,cl] z,^p„ (2.1) 

V2 ,j 
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where Gp is the Fermi constant and V stands for the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix. The operators Pj can be found in Eqs. (2), (73) and (101) of Ref. [9].^ The ones that 
are relevant for our present matching computation read 

Pi = {sLl,T-CL){cLrT%L), 
P2 = {sL7l.CL){cLYbL), 

q=u,d,s,c,b 

Ps = Z;'^{-SLa^''T%j,)G';^^, 

Pii = isLl,lu'ypT^CL){cLYYYT%L)-lQPi- (2.2) 
Their Wilson coefficients can be perturbatively expanded as follows 

^ C?^^) + a^C?^')+alC^('^ + alCf^'^ + 0{at), Q = c,t . (2.3) 

where 5s = as/(47r) = 5'^/(47r)^ and C"^*^"-* originate from ra-loop matching conditions. We 
neglect the 0{ae^) corrections to the r.h.s. of the above equation as well as additional operators 
that arise at higher orders in the electroweak interactions. 

The goal of the present paper is finding C^^^^ and Cg^^^ at the renormalization scale 
/io ~ {rrit or Mw)- As we shall see, it is convenient to consider different scales /iq for Q = c 
and Q = t. This is the reason why we refrain from applying unitarity of the CKM matrix 
throughout the paper. 

The renormalization constants Zij that enter Eq. ()2.H) are all known to sufficiently high 
orders from previous calculations [14,15]. The ones that are necessary here read (in the MS 
scheme with D = 4 — 2e) 

7 _ 116?r;2, 23848 i 13390 \ ~3 i (n(nA\ 7 _ 19 ?i;2 , 7249 , 5749 \ ~3 1 n(7iA\ 



81e s I 2187e2 ' 2187e / s ' ^V^s^' 27e s ' 1458e2 ' 5832e / s ' ^V^s^' 

^47 = -^^^l^O{al\ Z,, = -'^&'^+0{al), 

Zjy = (§ + 11)52 + 0(53), Z,s = 0, 

Zs7 = -ia,+ {B-m)^s+0{a!), Zss = 1 - f 5^ + (i| + H ) 5,^ + 0( 



Their overall signs correspond to the following sign convention inside the covariant derivative 
acting on a quark field ip: 



D,^ =[d^ + tgGlT- + leQ^AA ^. (2.5) 



^For simplicity, we set Vub to zero here, which makes irrelevant the operators P^" from Ref. [9]. The operators 
from that paper are denoted by Pj here. Our final results are insensitive to whether Vub vanishes or not. 
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For completeness, one should also mention the MS renormalization constant for the QCD gauge 
coupling in the five-flavour effective theory ((?bare = Zgg) 

Z, = l-i&s+{m-i)^! + <^i^!)- (2-6) 

Following Ref . [9] , we ignore the quark- mass and wave-function renormalization constants in the 
effective theory^ because their effects cancel anyway in the matching condition with analogous 
contributions on the full SM side. Only the top-quark contributions to these renormalization 
constants will be included in the SM counterterms (see Section H)). 

The coefficients Cl^""^ vanish for k = 1,2,11. At the tree-level, only C^^^^fio) = —1 is 
different from zero. All the ^^Vo) and ^^Vo) 

were found in Ref. [9] up to 0{e) and 

0{1), respectively. In particular, C2^^\fio) = and 



Ct'"{f,o) = -15 



w 





3 !i 






1 




~ 3 







-9x2 ^ ^g^, _ 4 _i ^ 21x2 - 42x - 4' 



6(x - 1)4 e 36(x - 1)3 ; 

/-45x2 + 38x + 28, 23x3 + 93x2 + 66x - 308 \ 

+ n 36(.-ir + — — 1+0(^,(2.10) 



where 



m2(/io) 



^ = (2-11) 
has been introduced. For later convenience we also define the variables 

w=(l-^^], z=^, y = ^^. (2.12) 

V m?(/io)/ ' ?«?(/^o) ' mtifio) 

In the following, the MS-renormalized top-quark mass mt{fJ^o) will often be denoted by just mj. 

For our present purpose, C^g^"* and C^g^"* are needed up to 0{e^) and 0{e), respectively. In 
practice, this implies a necessity or repeating the one- and two-loop matching computations for 
these coefficients from scratch. We shall describe this calculation together with the three-loop 
one in the following three sections. 

^although their non- vanishing vahies were relevant in the calculations of Zij H2.4(l and Zg H2.6(l 
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Figure 2: One-loop IPI diagrams for h 
background field gauge. 



S7 in the SM. There is no iy^0^7 coupling in the 



3 The unrenormalized SM amplitudes 

We have to consider all the one-, two- and three-loop one-particle-irreducible (IPI) diagrams 
contributing to the processes b —>■ s'j and b ^ sg. The one- loop b —>■ sj diagrams are shown in 
Fig.|^ Higher-order diagrams are found by adding internal gluons together with loop corrections 
on their propagators. 

We use the 't Hooft-Feynman version of the background field gauge for the electroweak 
interactions and QCD. Before performing the loop integration, the Feynman integrands are 
Taylor-expanded up to second order in the (off-shell) external momenta, and to the first or- 
der in the 6-quark mass. Thus, effectively, the only massive particles in our calculation are 
the top quark, the W boson and the charged pseudogoldstone scalar (j). The amputated IPI 
b ^ S'J Green function can be cast into the following form: 



AGf ePn 
' V2 (47r)2 



withPB= (l + 75)/2. 



y;{x) 



n>l 



Tie ~n— 1 



a, 



(3.1) 

(3.2) 
(3.3) 



n>l 



A. 



mS7 



e and At = where 7 is the Euler constant. The symbols 5*^ stand for dif- 

ferent Dirac structures that depend on the incoming 6-quark momentum p and on the outgoing 
photon momentum k 

[li^H^ 7m (P- ^), 7MP^ 7Mfc^ i)V^, l^p^, l^k^, 

nib^-ff,, mb'jf,^, rrib^-ff,, rrib-ff,^, M^7^) . . (3.4) 



5, 



The first two terms in the expansion of Y^^ ()3.2|1 are independent, but the third (three- loop) 
and higher terms do depend on x. 

By analogy, the b sg Green function reads 



AC nP I ^'^ ^'^ 



(3.5) 



i=i 
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with 

\nt ~n-l /^c{n) 



= E^r5r^Gf^(x), (3.6) 

n>l 

G]{x) = j:Arar'Gf-\x). (3.7) 



n>l 



As shown in Refs. [9, 16], only the following linear combinations of Y^^^^ and Cf^"'"' are 
sufficient for finding the coefficients C'j{^q) and C^{^q): 

c?!£.. ^ + (3.8) 

cZ^r. - Ig^-^ + gT- (3.9) 

The calculation of C^baL ^(^) requires supplementing Eqs. (57) and (58) of Ref. [9] 

by higher orders in e, which yields^ 

t(2) / -, (1 + 2e)jv-5(l + e)na+n3-3(l - e)l-n2(l ~ ^)l-n3 ^(1 + 2e)r(l - e) /o 

"™ -2=0 ^ ^ (^^_ 1)1(^2- 1)1(^3 -l)!(l-e) r(l + e) ^ ^ 

and 

= 2(r34r3i3{-^ + |'-i-'' + (i-2'-)i'>^'- + 2(i-'-)L..(i-i 



+2e(l 



2(i_e)(l-2e) 

1\ ..A 1 



Lig (1 - r) - Lis (l - -) - Li2 (l - -) Inr + e - InV + 0(e2)|, (3.11) 



for the generic two-loop integral 

r(2) ^ y^l) f " gl " g2 

n™ 7r4-2e r(l + e)2 7 (g2 _ ^2)„i (^2 _ „^2y„, [(^^ _ ^^)2]n3 ' J 

where r = m^/m^, = ni + ^2 + ^3 and (a)„ = r(a + n)/r(n). Otherwise, the calculation 
proceeds precisely as described in Section 5 of that paper. The unrenormalized one- and two- 
loop results read 

^c(l) _ 23 I 145e , 875£^ , 23£^7r^ , (o i o\ 

^7,haie — 36"^ 216 1296 432 ^ {O.IO) 

<L = i + t + W + ¥ + (3-14) 

r"'(2) _ 112 107 4147e , 56en^ , rOf'^^N /o i t-\ 

'^7,bare " 81e 243 1458 "^ 81 '^^\^)^ {■i.-LO) 

^c{2) _ 23 I 833 I 13429e , 23e7r^ , ^f{^ 

<^8,bare - 27i + 324 + ImT + + '-^V^ W'-^OJ 

ritC^) _ I f 3x''^-2x^ x'-l I 22x''-15:-i,i - + 159,x-46 

^7 harp — I J- "T ir, 



7,bare ^ 12 ^ V 4(x-l)4 e ^ 72(x--l):-' 

^All the other equations in Section 5.1 of Ref. [9] are vahd to ah orders in e. 
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-18a;3+150x^-157x+46 , 122a:='-933x^+975x-290^ 



72(x-l)'' 

-108x3+918x^-977a:+290 
432(x-l)4 



+ 



432(x-l)'^ J 
694x^-56193;^ +5937X-1750 



2592(x-l)3 



(3.17) 



24{a;-l)3 



^15x^-14x+8 I 13a:3+15a:^+186x-88 

24(a;-l)4 e 144(x-l)3 

35x3+273x^+1110x-680 



'^7,ha.ve 



I 2 / -81x2-130x+88 Ir, ^ _|_ 
+^ V 144{x-l)4 + 



864(x-l)3 



(3.18) 



-6x"t-46x^+28x^ _|_ 34x'^+101z»+402x^ -397X+76 A 



3{a:-l)5 



27{a;-l)4 



I -163:''-122x3+80x^-8x rr/ A , -333a;'* -2529x3 +688a:^ +778a:-224 x"-! 
9{x-l)4 n\^X,t)-\- 81(x-l)5 e 



-220x-*+12952x3-9882x^+2407x-397 
243{x-l)4 



+ e 



146x'*-4289x'='+2736x"'+14x-224 t ; / 1 
81{x-l)4 1 i 



'-^8,bare 



I -879x-*-50319x''+35810x^-5884x+428 ^„ ^ j_ -4381x-*+148252x3-89391x^+8797x-745 l , Ci 1 Q\ 

i 486(x-l)5 ^^^"^ 1458{x-l)4 J+^-^ltJ^ [6.1^) 

— 1 { ^ -L £2ll\ ( 17x3+31x^ x^-l _|_ 35x'* - 170^3 -447x^ -338x+56 A 
~ eV 6 ) \ 2(x-l)5 e 36(x-l)4 ^ 

. -4x'*+40x3+41x^+x o-/ \ _|_ -144x'*+4707x3+8887x^-122x-368 x"-l 
6(x-l)4 y-^-' ^> ^ 216{x-l)5 e 

641x4+184x3 +8001x2-220x-368i 



I -1367X-4 -9646x3 -76869x^ +7442x+2680 , 

1296(x-l)4 ^ 



I 2982x'*+30843x3+147437x^-7846x-6664 ]„ ™ _i_ 
1296{x-l)5 mx-^ 



216(x-l)'* 

22703x^-56674x3-934701x2 -46090X+59656 



Li2 1 - - 

X 



7776(x-l)4 



+ 0(e2), (3.20) 



where 



H{x,e) = Li2(l-i)+e 



Li,fl-a;) 



Lis (l - ^) + Li2 (l - \) Inx + iln^x] . (3.21) 



In addition to the bare coefficients, we shall also need those parts of Cf\^^^^ that originate from 
the mf,-dependent Dirac structure S'lo, as they play a separate role when gets renormalized. 
They read 



B7 



t{i) 



10 



1 I ( -3x2+2x x"-l I 5x2-3x \ , ( -2x2+x x"-l , llx2-5x 

12 A 6(x-l)3 e 12(x-l)2y' V4{x-1)3 e ^ 



24(x-l)2 



, 2 -12X-+5X 1 , 23X--9X , . 

+^ V 24(x-l)3 "^"^ + 48(x-l)2 ) + 



2 3x2 -9x 
(x-1)^ 



R = /^'(l) _ Z' 1 I e^TT^ ^ / X x^-1 I x2-3x ^ \ c ( 3x x^-1 I x2-7x ^ 
-D8 — LtiO - 12 A2(x-1)3 e 4(x-l)2y' ^ V4(x-1)3 e 8(x-l)2 ^ 



+^^(8(^1^^ + 



8(x-l) 



X — 15x 
16{x-l)2 



+ 0{€'). 



(3.22) 



(3.23) 



The renormalization of will not matter in the charm sector because 



c(l) 



10 



c(l) 
10 



0. 



Let us now turn to the main purpose of our paper, i.e. to the three-loop calculation. One 
of the 0(10'^) diagrams that we have calculated at this level is shown in Fig. ^ Obviously, 
when the virtual top quark is present in the open fermion line, we have to deal with three-loop 
vacuum integrals involving two mass scales, mt and My/- However, such double-scale integrals 



7 



are encountered in the charm-quark sector, too, when closed top-quark loops arise on the virtual 
gluon lines. 

At present, complete three-loop algorithms exist for vacuum integrals involving only a single 
mass scale. We have reduced our calculation to such integrals by performing expansions around 
the point rrit = Mw and for rrit ^ Mw- In the latter case, the method of asymptotic expansions 
of Feynman integrals has been applied [17]. At the physical point where Mi^/uit ~ 0.5, both 
expansions work reasonably well (see Section IHl). 

Two different approaches have been used for the calculation of the three-loop diagrams. 
The first one is based on a completely automated set-up where the diagrams are generated by 
QGRAF [18], further processed with q2e [19] and exp [20], and finally evaluated and expanded in 
e with the help of the package MATAD [13] written in Form [21]. MATAD is designed to compute 
single-scale vacuum integrals up to three loops. The individual packages work hand in hand, 
and thus no additional manipulation from outside is necessary. Moreover, all the auxiliary files, 
e.g. make-files to control the calculation or files to sum the individual diagrams, are generated 
automatically. 

The program exp is designed to automatically apply the rules of asymptotic expansions in the 
limit of large external momenta or masses. Thus, its output crucially depends on the limit we 
consider. For the expansion around = M\y, the asymptotic expansion reduces to the usual 
Taylor expansion in powers of w = (1 — M^/m^) and thus exp essentially rewrites the output 
of q2e to a format suitable for MATAD. However, for mt ^ Mw-, next to the Taylor expansion in 
z = M^r/mf, more diagrams expanded in various small quantities contribute according to the 
rules of asymptotic expansions. The package exp provides a proper input for MATAD which then 
performs the expansions up to the required depth, and computes the resulting scalar vacuum 
integrals. The mass scale of the latter is either given by rrit or M\y. 

An important element of the calculation are the so-called projection operations that pick only 
the two Dirac structures we need (see Eqs. (j3.8|) and (|3.9p ). and thus bypass the time-consuming 
tensor algebra. 

Using this method, we evaluated the expansions in z and w up to orders and w^, respec- 
tively. Furthermore, it was possible to compute the first few expansion terms for general gauge 
parameter ^, in order to check that it drops out in the sum of all bare three-loop diagrams. 
This imposes a strong check on the correctness of our results. 

In the second approach, MATAD was also used for three-loop scalar integrals involving a single 
mass scale. However, the diagrams were generated using FeynArts [22]. The remaining part of 
the calculation was performed with the help of self-written programs, largely overlapping with 
those used several years ago for the calculation of three-loop anomalous dimension matrices [14]. 
No projection operations were used, and all the Dirac structures (except for the ones quadratic 
in k) appeared in the results, which allowed for performing several consistency checks. This 
approach was obviously much slower, and was finally brought through thanks to the use of the 
Z-Box computer^ at the University of Ziirich. Only the expansion around mt = Mw (up to w^) 
was calculated using this method. 

Although our results for the three-loop diagrams are known in terms of expansions only, 

"^http : / /krone . physik . unizh . ch/^stadel/ zBox, 
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we are able to determine the exact x-dependence of their pole parts by using the matching 
equation discussed in Sectional Of course, we have verified that these pole parts have precisely 
the same expansions in z and w as found from the direct calculation up to and . 



Our results for C^bai-e C^baL ^^^e the following form: 



c(3) 
bare 



c\ 



a 



.c{3) 
8, bare 

<3) 

7, bare 

*(3) 

8, bare 



10798 1 
2187 e2 ^ 7 

4675 1 
1458 e2 ^ 7 

V\{x) + f\{x) 
V\{x) + fl{x) 



215 



224 



162 
11783 



In X 



648 



405 X 243 

169 46 , ■ 

1 ma; 

2160 X 81 



+ /s (^) 



where the pole parts in the top sector read 



-57xS+634x"^+1911x3-1044g^-4x 



9(x-l)6 



(in X + f In^ 



X 



+- 



380a:'^ -1099a:-^ -8521a:^ -4385a:^+5797x-812 
54{a;-l)5 

-560a;S+190a:'^+12410x^-6200a:^+496x 
27(x-l)5 



^Li2 1 - - 

X 



I -35586x5+223524x-*+1345261x3-604386x^-55425a:+20852 i 

1458{x-l)6 ■'"•^ 

I -325132x6+7070681x5 -72622435x"^+48723685x^-11218745x^+1543882x+864 



43740 x(x-l)5 
-199x"' -3018x3 -2535^2 -8x (\ . , e i„2 



12(x-l)6 



(in X + I In^ 



I 2054x5-11080x'*+52535x3+105505x^+26875x-3089 
360(x-l)5 



+ - 



+ 



-140x5+964x4-4813x3-3440x2-59x 



Li2 ( 1 - i 

X 



18(x-l)5 

-75843x"'-11835x-*-9946790x3-8078850x^+114225x+114293 i 

9720(x-l)6 ^^^-^ 

-5561837x''+4392955xS+397608280x"^+760910570x3-79703785x^-4603813a:+45630 

583200 x(x-l)5 



For the UV- and IR- finite functions /^(x), we write the expansions as follows: 

ln™x 



fcQ 
nm 



X" 



:-l)'™;z" ln"^z, 



(m < 3), 



1 



The values of and 6^*^ that we have found are listed in Appendix A. 



(3.24) 

(3.25) 

(3.26) 
(3.27) 



(3.28) 



(3.29) 

(3.30) 
(3.31) 
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4 The SM counterterms 



The renormalization scheme that we apply on the SM side is chosen in such a way that values 
of the renormalized as, the light-quark wave-functions and masses overlap with their MS coun- 
terparts in the five-flavour effective theory. Thus, as means af^jjio) throughout the paper. The 
one-loop renormalization constant of the QCD gauge coupling in the SM reads (cf. Eq. ()2.6|) ) 



Zf' = l + 'fl{--^^ tN, + 0{al). (4.1) 



5s / 23 1 

•••■ = ' + T(-ir + 3- 

Here, parametrizes the one-loop threshold correction that arises in the relation between af^ 
and af\ i.e. when the top quark is decoupled from as- A collection of explicit expressions for 
such corrections (also called "decoupling constants") up to three loops can be found in Ref. [23]. 
The value 

iV.= r(l + 6) (4.2) 

is found (exactly in e) from the requirement that the top-quark loop contribution to the off-shell 
background gluon propagator with momentum q is renormalized away, up to effects of order 
/m^ that match onto higher-dimensional operators in the effective theory. 
The same requirement applied to the light-quark propagators at two loops leads to the 
following expressions for the renormalization constants of their wave-functions and masses 

(^bare = 2"^^, "^bare = Zmm) 

AZ„ ^ ^SM_^e..tw. ^ fisX^(-^ + ^-|)+0(Ce). (4.4) 

In our calculation, the latter renormalization constant matters for the 6-quark only, because we 
include linear terms in m;,, while all the other light particles are treated as massless. 

The differences AZ^ and AZ^ are everything we need to know about the renormalization 
of the light-quark wave functions and masses. Since the wave-function renormalization matters 
for external fields only, the remaining parts of the considered renormalization constants cancel 
out in the matching equation, i.e. in the difference between the full SM and the effective theory 
off-shell amplitudes. It is worth noticing that since AZ^ and AZ^ arise at 0{a1) only, they 
had no effect on the two- loop matching computation in Ref. [9]. 

As far as the top-quark mass is concerned, we renormalize it in the MS scheme, at the scale 
/io, in the six-flavour QCD. The corresponding renormalization constant, when expressed in 
terms of fis = a^^-'(/io), takes the following form (exactly in e) 

4 / 74 97 S \ 

Z^,^l + a^Z^^l + + 0{al) = 1 - -5, + - _ - _iV,j al + 0{al). (4.5) 

Two more QCD renormalization constants need to be thought about in the context of our 
calculation. The first of them is the external gluon wave-function renormalization constant in 
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the b sg case. In the background field gauge, it just cancels with the renormalization of 
the gauge coupling in the vertex where the external gluon is emitted.^ The second one is the 
renormalization constant of the QCD gauge-fixing parameter ^. It plays no role either, because 
C?b!L and C^l2e are ^-independent.^ 

Last but not least, one needs to consider possible electroweak counterterms. Since we work 
at the leading order in the electroweak interactions, the only electroweak counterterms that 
may matter for us must have the sb fiavour content. Their dimensionality cannot exceed 4, 
and they must be invariant under the QCD and QED gauge transformations. These conditions 
leave out only two possible electroweak counterterms: sp b and sb. They originate from the 
fiavour-off-diagonal renormalization of the quark wave-functions and Yukawa matrices. Since 
we refrain from applying unitarity of the CKM matrix (but set Vub to zero), we write the 
corresponding electroweak counterterm Lagrangian as follows: 



V2 



0,sb 



(4.6) 



with the factors Ac and At that have been defined below Eq. ()3.3|1 . The renormalization 
constants and Z^^f^ are fixed by the requirement that the renormalized off-shell light- 
quark propagators with momentum q remain flavour-diagonal, up to effects of order /M"^ that 
match onto higher- dimensional operators in the effective theory. A simple one-loop calculation 
gives 



'2,sb 



'0,sb 



2 - 2e 



2-e 



0, 

rfe) 



X _ 2^ ~|~ 3^ 2 . ^ ^ > 3^ X 2 

IH z - 1) + e ^ — 

2 2{x-lf ^ ' 4(x-l) 

5/4x2-s + 2, -7x2-a;-2\ 
+e 1^7 7TT^lnx + — 7^ +C(e3) 



4(x- 1)2 



x-V 



x{x^ — x)T{e) 
(l-e)(x-l)- 



(4.7) 
(4.8) 



(4.9) 
(4.10) 



Higher-order (in a^) contributions to Z2sb and sfe are irrelevant to us, because the coun- 
terterms ()4.(jj) affect our calculation only when inserted into two-loop diagrams containing 
top-quark loops on the gluon lines. Otherwise, the loop integrals vanish in dimensional regu- 
larization after expanding them in external momenta, because all the propagator denominators 
are massless. As far as the tree-level diagrams are concerned, they give no contribution to the 
relevant structures 5*2 and Siq in Eq. (j3.4|) . 



'"^In the usual (non-background) 't Hooft-Feynman gauge, we would need to introduce, by analogy to Eqs. H4.3|l 
and A {Zg^/Z^ = a^N^ (-3/(4e2) + 5/(8e) - 89/48) + 0{a^,,e). 

^Contrary to the bare two-loop Wilson coefficients of the EOM-vanishing operators (Eq. (73) of Ref. [9]). 
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5 Matching 



We are now ready to write down the matching equation that follows from the requirement 
of equality of the effective theory and the full SM amputated IPI Green functions. The for- 
mer ones originate from tree-level diagrams only, because all the loop integrals with massless 
denominators vanish in dimensional regularization, after expanding them in external momenta. 

For the coefficients {i = 7, 8), the matching equation up to three loops takes the following 
form: 



n=l 



, 2(n-l) 



+S''^AIAZ^B, + \A\^A^Q al [z^.K, + Z^.R^) + 0{al). 



(5.1) 



Non-vanishing contributions on the l.h.s. arise for k = 1,2,4,7,8,11 that we have considered 
in Section |21 The effect of mb-renormalization is contained in the AZ^-Bj term, where Bi have 
been given in Eqs. (IT^ and (IT^ . 

The quantities t/"^ originate from the top-quark mass renormalization. Replacing in the 
bare results mj by Z^tf^t and Taylor-expanding in 5s, one finds tY = 0, 



.(2) 



.(3) 



2A^Z(i) (t— - f\ C'^^'> 



(5.2) 



2Z^D 



d 



2x' 



[1 - 4e)a; 



d_ 
dx 



2e' 



C 



t{i) 



,bare 



d_ 

dx 



+2ArzW[x^-2e]cl'^L. 



(5.3) 



The explicit factors of e in the above equation are due to the fact that At depends on rrit, too. 

The quantities Ki and Ri on the r.h.s. of Eq. (j5.ip originate from two-loop b ^ s'y and h ^ sg 
diagrams with insertions of the electroweak counterterm (|4.6|) and with closed top-quark loops 
on the gluon lines. We find 



^ , 88^ , Q^2^ 
405 ^ 6075 ^ )i 



Rr 



169 



10333 e 



2160 ' 64800 
1 I 193 8441 e 



40e ~'~ 1200 



36000 ^ 240 



(5.4) 



t(3) 
7 

t(3) 



It is interesting to notice that the mf,sb counterterm from Eq. ()4.fj|l is irrelevant for G 
(because R7 = 0{e^)) and for the charm sector (because Zg^^ = 0). Thus, it matters for C. 
only. 

At this point, all the ingredients of the r.h.s. of the Eq. ()5.H) have been explicitly specified. 
As far as the l.h.s. of this equation is concerned. Section |21 provides us with all the necessary 
renormalization constants and Wilson coefficients, except for and Cg . Thus, we can find 
Cj^^") and Cg^^"^ forn = 1, 2, 3 by solving our matching equation (j5.ip order- by-order in ctg. All 
the 1/e^ and 1/e poles cancel during this operation, as they should. The resulting finite Wilson 
coefficients are presented in the next section. 
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6 Results 

Our final results for the renormalized Wilson coefficients of the operators P^ and are as 
follows: 

- | + ^(il + |ln4)+e^(t + g + liln4 + lln^4)+0(e3), (6.2) 

r'c(2)/',, N _ _713 _ 4 1 , /^_7357 , STtt^ _ 820 Mo i HO l„2 \ , /'/^/',2x 

I/^OJ - 243 8T^"m^ + ^(^ T458 + ~8r ^ ^ ) ^ h \^-^) 

Cf\l^o) - -l + *ln4 + e(ii + ¥ + fgln4- + iln^^)+0(e^), (6.4) 
Cf\f,o) = C^fH/^0 = M^) + ^ln^ + filn^4- + O(e), (6.5) 



Cf\f^o) = CfH/.o = M^) + Sln^ + |iln^^ + 0(e), (6.6) 



5832 ' 486 



^ _ A I ^In '^O I 1„2 Mo , \ / 3x^^-2x2 , 22a»-153x^+159x-46 

<^7 l/^oj - (^i + e^n^ + yln + — J (^4(^_i)4 ^ + 72(x-i)3 

A _L ^ In /^o ^ f - 18a;^ +150x^ - 157a;+46 , 122a;^-933a:^+975a;-290 \ 

+e(^i + eln^j 72{x-i)i — + 432(^-1)^ ) 

I ^2 / -108x^+9180:^ -977a:+290 !„ ™ i 694a;3-5619a;^+5937a;-1750 \ , /n/'-r^^ ^fi 7^ 

"^^ V 432(x-l)4 mx-h 2592(1-1)3 J^'-^K^)^ K'^-') 

r^W(,, \ - (^ I rin 4 I ^0 I <'^'''^\ ( -3^^ ^''1 I 5a:'^-9a:^+30x-8 \ 
<-8 l/^oj - + + y m ^ + l^4(^_i)4 ^ + 24(^-1)3 j 

I ^ /i _i_ ^ l„ Mo \ / -15a;^-14j;+8 a^-l , 13a:3+15x^+186j;-88 \ 
+e (^i + eln^j I, 24(a;-l)4 — ^ 144(x-l)a ) 

I ^2 / -81x^-1303: +88 !„ „ i 35x3+273a;^+1110a;-680 ^ , (P. S^l 

V 144(a;-l)4 X + 864(x-l)3 )^^\^)^ 

-^t(2)/, \ -16x^-122j3+80x^-8x ^\ ( ^ , o^i^ Mn \ i -387a;'*-1413a;3+997a:^-65a:+4 x^-l 

W l/^oj - 9(^31)4 ii(x,ej (^i+^eln^j H 8i(r.-i)s ^ 

" UGx" -4289^3 +2736x^ +14J-224 _ 



_l_ 943:'^+186653;3-20682x^+9113x-2006 , ^ 
486(a;-l)4 ^ 



81(a;-l) 



I -1203j:^-433533:3+3703l3:^-1053l3;+1640 ]„ ™ i -6128a:'^+2528393:3 - 183912j;^ +43607 j;-7910 " 
486(x-l)5 niXi- 2916(a;-l)4 

_l_ ^In ifi / -720a;-^-3942a;3+1685a:^+713a;-220 !„ ™ i -346a:'^+44569a:3-40446a;^+13927a:-2800 A 
81(a;-l)5 mx-h 486(a;-l)4 j 

_l_ /'in I 3e i„2 Mo I e7r2\ ^ -6x-^-46a;3+28a:^ x^-1 , 34J;4+101a:3+402a;2-397x+76 ^ , n( (P. Oi\ 
"^V ^^^"^2^" 12 A 3(x-l)5 e ^ 27(x-l)4 lO-^; 

r'*(2)/',, \ _ -4a;-*+40a;3+41a;^+a; tr/^ ^\ /'i i o^ln ifi^ _i_ - 144a;-* +3177a;3 +3661a;^ +250a;-32 x^-l 

^8 KM - 6(^::T)4 -«l2;,ej (^i+^em^j H 2i6(^iap T' 



I 247j;^ - 1 18903:3 - 3 1779x^ +2966x- 1016 , 

1296(3:-1)4 ^ 



■t ■ 

64l3:'^+184a;3+8001g^-2203:-368 T ; Z' 1 1^ 
216(a;-l)4 ^^-^X} - x) 



_l_ 29823:4+22581x3+109751x^-1018j;-2968 !„ „ _i_ -18557a:'^-38590a;3-661839x^-100078j;+31096 
1296(1-1)5 mx-^ 7776(a;-l)4 
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, ^Yn-^ ( -'<'2a:'*+1971x3+31373;^+32j-100 ^ _|_ - 140a;^ -2692a;3 _ 13581a:^ +1301 j+208 N 
\ 54(x— 1)5 162(a;— 1)^ / 

_L /"in ^0 I 3e 1^2 i e^r^ \ /l7xf+31x^x^ , 35a;^ - 170a;3 -447x^ -338 j+56 ^ , n(^'^\ ( P, 1 fl^l 

V ^^^"^ 2 ^^^"^ 12 A 2(a;-l)5 e ^ 36(a;-l)4 ) ^ {O.IU) 

Cf\no) = Cf^ {no = m,) + In ^ [ -592.^-22.^+12814.3-6376.^+512. ^i, (1 _ l) 

I -26838.5+25938.^+627367.3-331956.^+16989.-460 i 

729{.-l)6 "^-^ 

I 34400.5 +276644.-^-2668324.3 +1694437.^ -323354.+53077 

2187(.-1)5 



_|_ £0. [" -63.5+532.^+2089.3-1118.^ x 
ml [ 9(.— 1)^ 

I 1186.5-2705.''-24791.3-i6099.^+19229.-2740 

162{.-1)5 



+ 0(e), (6.11) 



n^i^)(n\ — n^i^hu — rrj \ I In ^9 [ -148.5+1052.4-4811.3-3520.^-61. 7- /-I l\ 
<^8 i/^0 j - 1^8 i/^0 - "7,t j + m ^ 1^ 18(.-1)6 - 

I -15984.5+152379.'*-1358060.3-i201653.^-74190.+9188 i 

1944(.-1)6 -^"-^ 

I 109669.5 - 11 12675.4 +6239377.3+8967623.^ +768722.-42796 

11664(.-1)5 



mi 



2 ... 



I 1295.5-7009.4+29495.3+64513.^+17458.-2072 
216(.-1)5 



12(.-l)e 

VO(e). (6.12) 
As far as the three-loop quantities C^'^'^\no = Mw), C|*^^'*(/xo = Mw), C^^^\iJ,o = nit) and 

t(3) 

Cg if^o = ^t) are concerned, the matching calculation described in the previous sections gives 
us expressions for their expansions at x — > 1 and x — > oo. Denoting, as before, z — 1/x and 
w — 1 — z, we find 

Cj^^\fXo = Mw) ^ 1.525 - 0.1165^ + 0.01975^ In z + 0.06283z^ + 0.005349^^ \nz 

+0.01005^^ In^ z - 0.04202z^ + 0.01535^^ In z - 0.00329^ In^ z 
+0.002372^^ - 0.0007910^Mnz + O(z^), (6.13) 
C^^^) {fXQ = Mw) ^ 1.432 + 0.06709W + 0.01257w^ + 0.004710w^ + 0.002373m;^ 

+0.001406w^ +0.0009216^^ +0.00064730^/;^ +0.0004779^/;^ +0{w^), (6.14) 
Cs^^\fXQ^Mw) ^ -1.870 + 0.1010^ -0.1218^ In ^ + 0.1045^^ - 0.03748^Mn^ 

+0.01151z^ In^ z - 0.01023^^ + 0.004342^^ In ^ + 0.0003031^^ In^ z 
-0.001537^^ + 0.0007532^^ In ^ + O(z^), (6.15) 
Cl^^\lio^Mw) ^ -1.676 - 0.1179W - 0.02926w^ - 0.01297w^ - 0.007296«;^ 

-0.004672^^ - 0.003248^;^ - 0.002389«;^ - 0.001831w;^ + 0{w^), (6.16) 
C^^^) {f^Q = rrit) ~ 12.06 + 12.93^ + 3.013^ In ^ + 96.71z^ + 52.73z^ In z + 147.9z^ 

+187.7z^lnz- 144.9^^ + 236.1^Mn^ + e)(z^), (6.17) 
C^^^) {f^Q = rnt) ~ 11.74 + 0.3642^ + 0.1155«;^ - 0.003145^;^ - 0.03263«;^ - 0.03528«;^ 

-0.03076«;^ - 0.02504w^ - 0.01985«;^ + 0(w^), (6.18) 
Cf'\nQ^mt) ~ -0.8954 -7.043Z- 98.34^2 -46.21^Mnz-127.1z^ 
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Figure 3: The coefficients C7 ^"^(/io) as functions of y = / mt{fXo) . The (blue) dot-dashed lines 
correspond to their expansions in y up to y^. The (red) dashed lines describe the expansions in 
(1 — y^) up to (1 — y"^)^. The (black) solid lines in the one- and two-loop cases correspond to the 
known exact expressions. The (yellow) vertical strips indicate the experimental range for y. 

-181.6^=^ \nz + 535.8;z^ - 76.76z^ \nz + C(z^), (6.19) 
Cl^^\^Q = mt) ~ -0.6141 - 0.8975W - 0.03492^2 + 0.06791ti;^ + 0.07966u;^ 

+0.07226u;^ + 0.06132u76 + 0.05096w^ + 0.04216u7^ + 0{w^). (6.20) 

While only numerical values of the expansion coefficients have been given above, their exact 
values can easily be found from similar expansions for the unrenormalized three-loop results 
(Appendix A) and from the formulae of Sections I2HS1 

In Figs. El and m the top-mass dependent coefficients Cl^"\no = nit) and Cj^*-^''(/io = My/) 
for z = 7, 8 are plotted as functions oiy = My/ /mt{^Q). The different choice of renormalization 
scales in the top and charm sectors allows us to avoid logarithmic divergences at large and, 
consequently, achieve better control over the behaviour of the expansions. This is the main 
reason why has been normalized to My/ in the charm sector and to in the top sector, in 
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Figure 4: Same as Fig. El but for C8^^"^(/io 



all our intermediate and final expressions.^ 

The variable y changes from to 1, i.e. both starting points of our expansions are present 
in the figures. Note the relatively narrow ranges of the coefficient values on the vertical axes. 
The large rrit expansions (up to y®) are depicted by the dot-dashed lines, while the expansions 
around = M-w (up to (1 — i/^)^) are given by the dashed ones. In the one- and two-loop 
cases, solid curves show the exact results. The vertical strips mark the experimental values for 
y that we take (0.488 ± 0.015) for /io = rrit, and (0.461 ± 0.015) for /iq = Mw 

Comparing the three curves in the one- and two-loop cases (the two upper plots in both 
figures), one can conclude that a combination of the two expansions at hand gives a good 
determination of the studied coefficients in the whole considered range of y. However, the 
expansion starting from y = 1 works somewhat better for the physical values of rrit and Mw- 
Most probably, including more terms in the the large rrit expansion could improve its behaviour 
around y = 0.5. 

^Apart from that, many of the top-sector expressions would be significantly longer if /ip was normalized to 
M\Y there. 



16 



Although we do not know the exact curves in the three-loop case, the same pattern seems 
to repeat. In fact, the charm-sector expansions perfectly overlap in the physical region. In the 
top sector, one can (conservatively) conclude that 

C?^)(/io = mt) = 12.05 ±0.05, (6.21) 
Cf\fio = mt) = -1.2 ±0.1, (6.22) 

which is perfectly accurate for any phenomenological application. Let us note that a change 
of Cj (yUo = ''^i) from 12 to 13 would affect the b ^ S'j decay width by only 0.02%, while a 
similar variation of Cg^^\j2Q = rrit) would cause even a smaller effect. 

For the three-loop charm-sector coefficients, the uncertainty from the expansions is smaller 
than the one from the experimental error in Thus, one can safely use Eqs. fl6.1Hj) - ()6.16j] 
as they stand, without any additional uncertainty. Accurate values in the range 0.4 < y < 0.6 
can also be found from the following fits: 

/ m \ 0-0338 

Cf\fio = Mw) = 1-458 , (6.23) 

, ^ / rn \ 0-0598 

Cf\f^o = Mw) = -1.718 . (6.24) 

It is instructive to study the behaviour of the three-loop top-sector coefficients in a plot 
where subsequent terms of our expansions are successively taken into account. This is shown 
in Fig. El The quality of the two expansions in various regions of y is transparent there. 




Figure 5: The three-loop top-sector coefficients. The solid lines represent the highest orders we 
know (as in Figs. Q] and IH). The dashed and dot-dashed lines show the lower orders. 
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7 Conclusions 



The three-loop matching conditions found in the present paper complete the first out of three 
steps (matching, mixing and matrix elements) that are necessary for finding the NNLO QCD 
corrections to -B ^ -^s7- The effect of the NNLO matching alone is scheme- and scale- 
dependent. In the MS scheme with Mw < hq < nit, if stays within 2% of the decay width, 
i.e. it is significantly smaller than the total higher-order perturbative uncertainty that was 
estimated in Ref. [6]. This uncertainty is expected to get significantly suppressed in the near 
future, after the remaining two steps of the NNLO calculation are performed. 

The methods that we have applied in the present work are, in principle, applicable to any 
three-loop matching computation involving several different mass scales. A detailed description 
that we have presented for each step of our procedure can serve as a guideline for treating similar 
problems in various domains of particle phenomenology. 
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Appendix A: Three-loop expansion terms 

In this appendix, we present our results for the coefficients and b^^ from Eqs. (j3.3U|l 
and (|3.31|) up to n = 4 and n = 8, respectively. They are given in terms of the following 
symbols (see also Eq. (16) of Ref. [13]): 




16 4 16 



45 TrVSin^S 
T 8 



3571^ 9 




C2 + Ca + 6 V3CI2 - - 6 Vsim Lis 



i7r/6 \ 



(A.l) 
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where C^ix) — Im [Li2 (e'^)]. 



The expansion coefficients that we have found read 

I 70C(3) I 158752 _|_ 437r^ , 41634l7r^ , 46.D3 _ 92B4 , 820640533 

243 280 405 612360 81 81 9185400 

, 307721C(3) _ &7TI _ 28432752 1 969597r-* _ 5388025l7r^ , 67;^ 

324 18 840 116640 816480 27 
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94143997C(3) _ 67637r^ _ 113487750073 

1244160 15120 979776000 ' 

91942073C(3) _ 24l7r2 _ 10092619036343 

860160 1512 
137418234607C(3) _ Stt^ 

891813888 135 
7490373009073C(3) _ 5937r^ 

35672555520 30240 23597395476480 ' 

467301421361C(3) _ 13l7r^ _ 6438103242654889429 

1698693120 47520 19467851268096000 ' 



76814438400 ' 
13709395882765691 
73741860864000 ' 
5960644239577417 



329068267226885C(3) , 11397r^ 

941755465728 249480 

1164930029277053C(3) , 3537r^ 

2690729902080 46332 



1798715398515135307759 

4282927278981120000 ' 
445971686554467633047 
857084922534297600 ' 



32688338029429333C (3) , 562937r^ _ 12555120069922446322011879253 
62185757736960 



6486480 



19873056681818251591680000 



(A.3) 



22301C(3) I 2214952 1 177r'' , 1706597r^ , I3D3 _ 13.B4 1 1189623529 
648 224 2160 25515 216 108 61236000 ' 
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147193C(3) I 1075273S2 _ 17337r^ , 317238l7r^ _ 56.D3 i II2B4 1 409119733 

432 2240 2430 544320 27 27 1360800 

581 



581 
5400 ' 





80 ' 

304559C(3) I 49rf 342275952 _ 92897r^ , 2719077r^ _ 495| _ 349D3 1 349.B4 1 487461187 

216 2 1920 4320 31104 3 24 12 155520 ' 

I 1323^2 _ 381l7r^ _ 8187697 

4 432 64800 ' 



63517 
2160 
199 
" 72 ' 



8128418C(3) I 3891425Tf 28450834752 _ 2616726177^ , 11597258571^ _ 38914255| 
2187 17496 36288 12597120 ~'~ 5878656 26244 

6O97-D3 I 6O97-B4 I 2 
108 "t" 54 



I 389142552 _ 44209l7r^ _ 1407902803 
~'~ 1296 3888 583200 

10526363 

19440 
117 
2 ' 



26244 

197-B4 I 255386869021 
54 14696640 



1201430399C(3) , 72196517Tf 28165051597 52 , 835449797r-^ , 339720929337r2 _ 721965175| 
139968 69984 933120 10077696 117573120 104976 

34039 J>3 I 34039-B4 4. 206714107565 
216 108 3359232 



I 7219651752 _ 92809857r^ _ 4347341779 
5184 15552 291600 

20676181 

6480 
1379 
4 ' 



3426427C(3) , 40077r^ , 383324521 

207360 2880 4665600 ' 

377401C(3) _ 90l7r^ _ 683934529 

12960 10080 16329600 ' 

2r)702o:-i()i(;(3) , 24771'' 5133539931187 



(A.4) 



3870720 ' 20160 64012032000 ' 
148678249549C(3) , 257r^ _ 590349164605337 



1486356480 ' 1008 4916124057600 
48644809387C(3) , 8697r^ _ 33797460512670161 



339738624 ' 40320 196644962304000 ' 
3067923823757C(3) _^ 36497r^ _ 6031703724292608407 



15854469120 ' 221760 25957135024128000 ' 
787852603809259C(3) _^ 161897r^ _ 9469141155025123085807 



3139184885760 ' 1330560 31408133379194880000 ' 

113217484992547C(3) , 1737r^ _ 1857857271774679586364179 

358763986944 19305 4899668807154401280000 ' 

96339777793582171C(3) , 114857r^ 948680668305509273934263213 ( A 

248743030947840 1729728 2038262223776230932480000 ' V^-^^ 

5032C(3) I 296837r^ _ 517861 
243 4374 39366 ' 

I 407 
"""729 ' 
112 
243 ' 
112 
1215 ' 

~405 ' 
I 197r2 _|_ 4960261 



5670 ' 166698000 



20 



*21 
"22 

cP^ - 

"30 — 

"31 — 
"32 — 
"40 

"41 



2123 
396900 ' 

19 
1890 ' 
47r2 



3645 
5876 
"382725 ' 
4 

"1215 ' 
3075421 
1296672300 
74 
"93555 ' 



752008 
24111675 



(A.6) 



hi' 
hi' 

"3 

1,7c 



hi' 
hi' 
hl'^ 



7048C(3) 

243 
448C(3) 

27 
224C(3) 

9 

896C(3) 

27 
1120C(3) 

27 
448C(3) 

9 

1568C(3) 

27 
1792C(3) 

27 
224C(3) 



69S2 I 1039257r^ _ 6854384 
14 15309 3444525 ' 

I 65352 _ 297r^ , 670214 
35 8505 42525 ' 

I IOIO9S2 _|_ _|_ 16246 

210 ~'~ 17010 ~'~ 945 ' 

I 3676S2 I 47r^ _|_ 7423 

45 ""^ 3645 ~'~ 405 ' 
_|_ 316652 I 92539 
27 4860 ' 

_l_ 6251252 I 352919 
405 18225 ' 

_l_ 520052 I 235619 
27 12150 ' 

I 197088452 I 1468651 
' 8505 ' 76545 ' 
I 1385980952 1 6903803 

~'~ 51030 ~'~ 367416 ' 



(A.7) 



"11 — 

"20 ~ 

"21 — 

"22 — 

"30 — 

"31 ~ 

"32 — 

"40 — 

"41 — 



"00 — 

"01 ~ 

8c _ 

"02 ~ 

"10 ~ 



1627C(3) _^ 255837r2 , 9148337 



162 
1049 



5832 



+ 



104976 ' 



486 ' 

23 
81 ' 
9707 



32400 ' 
601 



+ 



2160 

297r2 
"7560 

2479 
'66150 ' 

29 
"2520 ' 

llTT^ 

" 108864 
496373 
114307200 ' 
11 



462793 
6945750 ' 



1616438903 
144027072000 ' 



36288 ' 
13282901 
"8644482000 ' 
1879 



2494800 ' 



(A.8) 



l8c 
'^O 

U8c 



2833C(3) _ 2987152 
162 2240 

131C(3) , 59752 
9 2240 



71696697r'' 
1632960 



7342571207 
73483200 ' 



14477r2 _|_ 17694343 



181440 



907200 
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78c _ 65C(3) I 50339S2 i 75l7r^ , 1260331 

"2 3 ~'~ 2240 ~'~ 181440 ~'~ 60480 ' 

78c _ 259C(3) I 64417S2 _ Utt^ i 4083773 

"3 9 ~'~ 1344 108864 181440 ' 

t8c _ 323C(3) I 10871S2 _i_ 307651 

4 9 ~'~ 144 12960 ' 

_ _AQ/-(Q\ I 4522352 i 4783799 

"5 - "i^^y^)^ 432 + 194400 ' 

l8c _ 451C(3) I 243IS2 _|_ 816431 

'^6 9 ~'~ 18 ~'~ 32400 ' 

78c _ 515C(3) _|_ 377316752 1 10434863 

"7 9 ~'~ 22680 ~'~ 408240 ' 

t8c _ 193C(3) I 10801258952 1 252135383 /a q\ 

"8 3 ~'~ 544320 9797760 ' l^-^J 
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